Nuclear recoil effect on the magnetic-dipole decay rates of atomic levels 
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Abstract 

The effect of finite nuclear mass on the magnetic-dipole transition probabilities between fine-structure 
levels of the same term is investigated. Based on a rigorous QED approach a nonrelativistic formula for the 
recoil correction to first order in m e /M is derived. Numerical results for transitions of experimental interest 
are presented. 
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Investigations of transition rates and transition energies provide direct access to most impor- 
tant spectroscopic properties of atoms and ions. In addition to the detailed understanding of the 
atomic structure, studies of transition rates in multicharged ions are of great relevance in plasma 
diagnostics and astrophysics [1]. The relative intensities of electric-dipole (El) forbidden transi- 
tions are frequently employed as a sensitive tool for plasma density diagnostics and coronal lines 
analysis. With the help of novel devices such as ion traps and Electron-Beam Ion-Traps (EBITs), 
new experimental high-precision data for the magnetic-dipole (Ml) decays have become available 
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The Ml -transition rate between the fine- structure levels (ls 2 2s 2 2p) 2 P ; 
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in B-like Ar 13+ ion associated with the 4412.4 A coronal line has been measured with an accuracy 
better than one part per thousand [|4j,|5|]. The elaborative theoretical calculation of this decay rate, 
including relativistic, interelectronic-interaction, and quantum-electrodynamic (QED) corrections, 
has revealed a discrepancy between theory and experiment QVL ISD - The total theoretical uncertainty 
is mainly set by the uncalculated recoil correction together with the experimental error of the tran- 
sition energy. In this work we present the evaluation of the recoil correction to the Ml -transition 
rates between fine- structure levels of the same term. 

The finite nuclear mass effect on the photon emission was first taken into consideration for 
El-decays in Ref. B9Q. The relativistic -recoil correction of the order (Zm e /M)(aZ) 2 (m e and M 
are the electron and nucleus masses, respectively) for these decays was calculated in work Jioll . 
The recoil effect on the two-El -photon transition 2s — Is and the transitions between hyperfine- 



structure levels was investigated in Refs. [ 



in light atoms were presented in Refs. [|13L 



1211. Some developments for the forbidden transitions 



14fl . However, the recoil effect on the Ml -decay 



between fine- structure levels has not been considered before. The derivation of the nonrelativistic 
formula for the recoil correction to the corresponding Ml -transition rates within a rigorous QED 
description is the main result of the present paper. 

Our consideration starts from the nonrelativistic Hamiltonian for N electrons and for the nu- 
cleus (atom) interacting with the second-quantized radiation field. Utilizing the transverse gauge 
and the Schrodinger representation it can be written as (in units h = c=l,e = — |e| < 0) 
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where the indices "e" and "n" designate the action on the electron and nucleus variables, re- 
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spectively, A(r) is the (transversal) vector potential of the quantized electric (£ t = —d t A) and 
magnetic (H = V x A) fields, and the operator V(r[ e ^ — r^) describes the Coulomb interaction 
between the i-th electron and the nucleus. The last term in Eq. (QQ) corresponds to the interaction 
of the electron spin s^ e ^ with the quantized radiation field, where /i — \e\h/2mc denotes the Bohr 
magneton. Here we assume that the nucleus possesses zero spin. It is most suitable to introduce 
relative and center-of-mass R variables for the electron-nucleus subsystem 



(2) 



R = [ Mr (n) + m e V r (e) | . (3) 

M + Nm e I i ) 

For the corresponding conjugate momenta and P one derives 

» = -ira («•'"' + (4) 

P = p<"> + £ vf • (5) 
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Performing corresponding substitutions and exploiting the condition Nm e <C M , the Hamiltonian 
(0Q) can be expanded up to zeroth order in a and to first order in m e /M. It may be written as the 
sum of an unperturbed (H ) and an interaction (H 7 ) Hamiltonians 

H = H Q + H 1 , (6) 

where H is defined as 

^i£P? + £Hr,) + i£^ + 40>) ! 

i i i^k i 

+ 2^M P2 + lJ d M^(x)+^ 2 (x)] (7) 
and H 1 describes the interaction with the quantized radiation field A(r), 

# 7 = -— p« • A ( R + r *) + s * (e) • n ( R + r *) 
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-^E p - A ( R +^) + § p - A ( R )-§E^- A ( R )- ( g ) 

According to the basic principles of quantum electrodynamics fislj . the S-matrix element corre- 
sponding to the transition of an atom from the state A to B accompanied by the emission of the 



photon with wave vector kf and polarization e/ is given by 

/oo r r , p2 P 2 \ 

dt / dR / dn . . . d riV exp (i^j|t + + ik° f t - i^t - iE A t) 

x^(R) ^(n, ...,r N ) (k f , e f \H 7 \0) * A (R) $ A (n, . . . , r w ) . 
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Here the wave function 



(10) 



describes the free translational motion of the atomic center-of-mass (it is supposed that the normal- 
ization volume equals to unity). The iV-electron wave function &a{ y i, ■ ■ ■ > r/v) = ( r i> • • • , r^y |^4) 
obeys the following equation of motion 
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i?4 and £"5 are the energies of initial and final many-electron states, respectively, m r = 
m e M/(m e + M) is the reduced mass. To evaluate the matrix element (kf, ef\H^\0) in Eq. © 
we employ the commutation relations for the creation and annihilation operators. Then in the 
spectral expansion of the quantized electromagnetic field only the term with wave vector kf and 
polarization e/ survives. In such a way, after the integration over t and R variables, the S-matrix 
element takes the form 

p2 



S~ff,B\A 

where 



-i (2tt) 4 5 3 (P B + k f -P A )5^ + E B + k° f -^- E A ) (B\V 7 \A) , (12) 
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A/ is the photon wave function, 



A f (r) 



Efexp(ihf ■ r) 



(14) 
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and Hf(r) = VxAj(r). One can notice that the last term in curly brackets of Eq. (fT3l) represents 
the leading recoil contribution of the order Zm e /M to the El-decay. In view of Eq. (fl~2l) the 
transition rate per unit of time can be written as 



dW 



S[^k + E B + k° 
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Taking into account the identity 



with/(x r , 
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and assuming = 0, one can easily find 
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where we keep only physically meaningful root and neglect the contributions of the order (Ea — 
E B )/M and higher. Expanding further the photon wave function in a multipole series (see, e.g., 



Refs 



.fl 



16(]) we keep only the Ml -component. The first three terms in braces of Eq. (1131) will 



contribute to the Ml -decay. Finally, the Ml -transition probability between the many-electron 
states \A) and \B) accounting for the recoil correction to first order in m e /M is given by the 
following nonrelativistic formula in terms of a reduced matrix element 



W = - 



4 u 3 f4 
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where the integrations over the photon energy and angles have been carried out, J a is the angular 
momentum number of the initial state, L and S are the orbital and spin angular momentum opera- 
tors, respectively, and uj = kj = Ea — E B characterizes the transition energy. The wave functions 
| A) and \B) are determined by Eq. (fTTI) . However, since we restrict our consideration to the first 
order in m e /M, the recoil corrections in Eq. CCD) can be omitted. The Ml-transition operator 
entering Eq. (TTSl coincides with the corresponding recoil corrected magnetic moment operator 



i,k 



(19) 



which was derived in Ref. fl 1 VTI for the evaluation of g factor. The calculation of the recoil correc 



tion to the bound-electron g factor to all orders in aZ can be found in Refs. tl8| 
Finally, the recoil correction to the Ml-transition probability takes the form 
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where \A) and \B) are linear combinations of the Slater determinants constructed in terms of 
one-electron Schrodinger wave functions and being eigenstates of the operators J 2 , J 2 , L 2 , and 
S 2 . This formula holds for all Ml -transitions between the fine-structure levels of the same term 
with arbitrary number of electrons. An additional recoil correction to Eq. (1201 arises also from 
the transition energy. However, we do not consider this contribution here, since we will further 
employ the experimental value for the energy of the emitted photon. 

In what follows we present some explicit results for highly charged ions restricting to the ap- 
proximation of noninteracting electrons. Let us consider first the Ml -transition between the levels 
(ls 2 2s 2 2p) 2 P3/2 — 2 - Pi 12 , m boronlike ions. To perform the angular integration we employ formu- 
las presented in Refs. 112(1 12111. Accordingly, for the second reduced matrix element in Eq. (1201 , 
we obtain 

((l S 2 2 S 2 2p) 2 P 1/2 ||^[r J xp fc ]|| 2 P 3 / 2 ) = -^=--^= (2D 

i,k V V c=ls,2s 

with 

n 2p , c = / dr l rfR 2p {r 1 )R c {r 1 ) / dr 2 r\ R 2p (r 2 ) — R c (r 2 ) , (22) 

JO JO dr 2 

and R c is the nonrelativistic Schrodinger wave function of an electron in state c. The radial inte- 
grations in Eq. (|22l) are performed analytically with the following result for the recoil correction 

AH- = ^^(l-^). (23) 

For the experimentally most interesting case of 40 Ar 13+ ion |4j,|5|] it yields AW Tec = 0.0017 s _1 , 
utilizing the experimental value for the transition energy u = 22656.22 cm -1 [22]. The theoretical 
value for the corresponding transition energy obtained, recently, by means of an ab initio QED 
approach |23|] is in perfect agreement with the experimental one. For this ion, the uncertainty of 
the recoil correction to the transition rate due to correlation effects is assumed to be about 20% 
by the relative estimation as a ratio n c /Z, where n c is the number of the core electrons. Adding 
the recoil correction to the decay rate calculated in Ref. Q8|] , we obtain the total values for the 
transition probability Wtotai = 104.848(1) s _1 and the corresponding lifetime r total = 9.5376(1) 
ms. The comparison between our theoretical result with the most accurate experimental value 
T xp = 9.573(4) (lg 2 ) ms (stat)(syst) Q3U reveals a remaining discrepancy of about 5.5a _ , where 
o~ is defined as combination of the statistical and negative systematic uncertainties in quadrature. 
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Furthermore, we consider the Ml -transition between the levels (ls 2 2s 2 2p 6 3s 2 3p) 2 Pzi2 — 2 P\/2 
in aluminiumlike highly charged ions. The angular integration of the corresponding reduced ma- 
trix element yields 

((l S 2 2 S 2 2/3 S 2 3p) 2 P 1/2 ||^[r i xp fc ]|| 2 P 3 / 2 ) = --^-^ ]T TZ 3PiC . (24) 

i,fc V V c=ls,2s,3s 

After the radial integration is performed the recoil correction results in the following expression 

9 Po M V 2 10 5 11 J 

For the case of 56 Fe 13+ ion, where the most accurate experimental result exists [60, we obtain 
Aiy rcc = 0.00057 s -1 . We assign a 50% uncertainty for this value due to the recoil correction 
beyond the one-electron approximation. The value of the corresponding transition energy uj = 
18852 cm -1 is taken from Ref. [24]. Finally, the total theoretical values for the decay rate and 
the lifetime are Wtotai = 60.44 s _1 and r to tai = 16.545 ms, respectively. These values include 
the interelectronic-interaction correction with the frequency-dependent term (see, for the details 
Refs. Ji |25 



]), the QED contribution, and the recoil correction. The details of this calculation 
will be presented elsewhere. The most accurate experimental value r cxp = 16.726+o;o 2 q ms Q6U 
significantly deviates from our theoretical result. 

Summarizing this work, the nonrelativistic formula accounting for the recoil correction in first 
order in m e /M to the Ml -transition rates between the fine- structure levels has been derived within 
a rigorous QED approach. Explicit results for the B-like and Al-like highly charged ions have 
been obtained within the one-electron approximation. The values for the recoil correction together 
with the total values of the transition probabilities have been obtained for the case of B-like Ar 13+ 
and Al-like Fe 13+ ions. Comparison of our total results with the most accurate experimental data 
reveals remaining discrepancies, which can not be explained by the recoil correction. However, 
we expect that further investigations of the transition rates in ions with several electrons could 
elucidate the origin of this disagreement. 
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